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We update our recent calculation of quarkonium Euclidean correlators at finite temperatures 
in a potential model by including the effect of zero modes in the lattice spectral functions. These 
contributions cure most of the previously observed discrepancies with lattice calculations, supporting 
the use of potential models at finite temperature as an important tool to complement lattice studies. 
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Quarkonia at finite temperature are an important tool 
for the study of quark-gluon plasma formation in heavy 
ion collisions (see, e. g., Ref. *lj). Many efforts have been 
devoted to determine the dissociation temperatures of 
QQ states in the deconfined medium, using either lattice 
calculations of quarkonium spectral functions @, OJ 0, d, 
IE S Hi or non-relativistic calculations based upon some 
effective (screened) potential @, HE HO, El M, QUE Eft 

Lattice studies are directly based on quantum chromo- 
dynamics and should provide, in principle, a definite an- 
swer to the problem. However, in lattice studies the spec- 
tral functions have to be extracted — using rather lim- 
ited sets of data — from the Euclidean (imaginary time) 
correlators, which are directly measured on the lattice. 
This fact, together with the intrinsic technical difficul- 
ties of lattice calculations, somehow limits the reliability 
of the results obtained so far, and also their scope, which 
in fact is essentially limited to the mass of the ground 
state in each QQ channel. 

Potential models, on the other hand, provide a sim- 
ple and intuitive framework for the study of quarkonium 
properties at finite temperature, allowing one to calcu- 
late quantities that are beyond the present possibilities 
for lattice studies. The main problem of the latter ap- 
proach is the determination of the effective potential: al- 
though at zero temperature the use of effective potentials 
and their connection to the underlying field theor y is well 
established, at finite T the issue is still open (TtI fl8l . 

Calculations of the cc and bb dissociation tempera- 
tures, using different potential models based upon the 
lattice free and internal energies, have found on the whole 
a reasonable a gre ement with the results from the lattice 
studies [Til. [l2l Il3l. [l5j . On the other hand, calculations of 
Euclidean correlators using a variety of potential models 
were not able to reproduce the temperature dependence 
of the lattice correlators 0, EE HH • 

In Ref. [l5| it was pointed out that a precise quanti- 
tative agreement with the lattice correlators should not 
be expected, mainly because lattice spectral functions are 
strongly affected in the continuum region by artifacts due 
to the finite lattice spacing. Actually, while some degree 
of qualitative agreement had been found for the S'-wave 



correlators, it was somehow disturbing the finding that 
for the P-wave correlators the temperature dependence 
of the potential model was even qualitatively different 
from the lattice one. In Ref. [lj| this different behav- 
ior has been shown to be due to a sizable amount of 
strength at low energy that was present in the P-wave 
lattice spectral functions, but not in the potential model 
estimate. This sizable amount of spectral strength - 
which is present below the QQ threshold even when no 
bound states are supported — strictly speaking goes be- 
yond a potential model description. Hence, at very low 
energy the effective potential approach has to be supple- 
mented by different mechanisms. 

Recently, in Ref. [2l| it has been shown that the lattice 
calculations of meson correlators at finite temperature 
contain a constant contribution, due to the presence of 
zero modes in the spectral functions. While the presence 
of a zero mode in the vector channel had already been 
discussed in the literature, in the P-wave channels it had 
generally been overlooked. 

In the following, we want to show that the inclusion 
of the zero mode contributions in the calculations of 
Ref. can reconcile potential model and lattice cor- 
relators. These contributions have also been considered 
in Ref. [161 ] and we shall comment below on the results 
obtained there. An alternative explanation of the poten- 
tial model vs lattice discrepancy has been proposed in 
Ref. [IH and, again, we defer at the end our comments. 

Here we briefly outline the model developed in 
Refs. [HI, [l5j], to which we refer the reader for all the 
details. The object we are interested in is the QQ Eu- 
clidean time correlator at finite temperature T, defined 
as the thermal expectation value of a hadronic current- 
current correlation function in Euclidean time r for a 
given mesonic channel H (see, e. g., Ref. [22j], Chap. 7): 

G H (r,T) = {j H (r)jU0)), (1) 

where j H = qT H q and T H = 1, 75, 7 M , 7^75. The 
four vertex operators Th correspond, respectively, to 
the scalar, pseudoscalar, vector and axial- vector mesonic 
channels, which in turn, at zero temperature, correspond 
to the Xco (Xbo), Vc (Vb), <//* (T) and \ci (Xbi) quarko- 
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nium states for the cc (bb) system, respectively. In some 
lattice studies only the spatial components in the vector 
and axial- vector channels are considered (Th — 7», 7i7s)- 
Moreover, we shall restrict ourselves, as in most lattice 
calculations, to the case of spatial momentum p = 0. 

The correlators of Eq. |T]) are related to the corre- 
sponding spectral functions through an integral trans- 
form, 



G h (t,T) 



du)a H (u,T)K(T,uj,T), 



which is regulated by the temperature kernel 

cosh[w(r- 1/2T)] 



sinh[w/2T] 



(2) 



(3) 



The temperature dependence of the correlators is usually 
studied by introducing the ratio between the correlation 
function Gh (t, T) at some temperature T and the so- 
called reconstructed correlator, 



GW(t, T, T r ) = / du <t h (uj, T r )K(r, 
Jo 



w,T), (4) 



calculated using the kernel at the temperature T and 
the spectral function at some reference temperature T r . 
This procedure should eliminate the trivial temperature 
dependence due to the kernel and differences from one 
in the ratio should then be ascribed to the temperature 
dependence of the spectral function. 

In the potential model of Rcf. [15], the spectral func- 
tion has been expressed as 



a H (uj, T) = J2 F^J(u - M n ) + 9(u - 



(5) 



where in the right hand side the sum over n runs over 
the bound states of mass M n and the last term repre- 
sents the continuum contribution, so being the contin- 
uum threshold; F H n and F H e are the couplings associ- 
ated to the discrete and continuum states, respectively, 
and they can be expressed [20|, |23| in terms of the wave 
function at the origin for the S-states (pseudoscalar and 
vector channels), 



f 2 ps = ^\m\ 2 



and Fi 



= ^l*(0)| 2 , 

Z7T 



(6a) 



and in terms of the first derivative of the wave function at 
the origin for the P-states (scalar and axial- vector chan- 
nels) , 



Fa 



9N C 
1-Kms 



| Ft (0)| 2 and 



9N C 



\R'(0)\ 2 , (6b) 



N c being the number of colors and m the quark mass. 
Since the integration over the energy in Eq. ^ can reach 
very high excitation energies, the wave functions have 
been renormalized to account for relativistic kinematical 
effects, that is the spectral functions have the correct 
asymptotic dependence (~ uj 2 ). 



The potential that we employ to generate the QQ wave 
functions [T3, EH is given by a linear combination [llj 
of the color-singlet free and internal energies obtained 
on the lattice from the Polyakov loop correlation func- 
tions [13, [H, H(| Ht| • The only partially free parameter 
one has in this approach is the bare heavy quark mass 
[30j : for instance, using the physical mass for the quark 
c, the J/tp dissociation temperature occurs below 1.5T C , 
whereas with slightly heavier quarks it can be brought 
above 1.5T C [l|. 

When addressing the Euclidean correlators, one should 
also account for the presence of zero mode contributions 
to the spectral functions at finite temperature. There is 
no zero mode in the pseudoscalar channel, whereas the 
vector channel has been discussed in Refs. [13, HH. Re- 
cently, is has been pointed out that these modes provide 
an important contribution also in the lattice scalar and 
axial- vector correlators [2l| . In the non- interacting case, 
they yield a term proportional to 5(u>) in the spectral 
function and a r-independent term in the Euclidean cor- 
relator: 



7 (o) 

-.(0) 



= Xh(T)lvS{uj) 



Gh = T X h(T). 



(7a) 
(7b) 



In the interacting theory, one may assume the delta func- 
tion to be smeared to a Lorentzian, with a width related 
to the heavy quark diffusion constant [28[ . 

Here, for simplicity, we employ the free expressions of 
Eq. ||7J), where all the susceptibilities \h can be expressed 
in terms of a scalar piece, 



dk m 2 drip 



J (27r) d uj£ duj k 



and a charge piece, 

X0 (T) - -AN C 



dk drip 
(2tt)3 8u k ' 



(8a) 



(8b) 



as 



pseudoscalar 


Xps 


= 


scalar 


Xs 




vector (7^) 


Xv = 


-Xs 


vector (7$) 


Xv = 


Xo - Xs 


axial (7^75) 


XA 


3xs 


axial(7i7 5 ) 


x • = 


Xo + 2%s 



(9) 



In Eq. flgj), uj k = Vfc 2 + m 2 and n F = l/[exp(w fc /T) + 1]. 

The zero mode contributions to the spectral func- 
tion arise from processes in which the meson current, 
rather than inducing transitions resulting in the creation- 
annihilation of a QQ pair, is absorbed by a heavy (anti- 
)quark of the thermal bath, which is in turn scattered 
with a sligthly modified momentum. The detailed bal- 
ance between this process and its inverse results in the 
factors driF/dujk in Eq. ([5]). 
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FIG. 1: Ratio of the charmonium S-wave Euclidean corre- 
lators at various temperatures to the reconstructed one at 
T — 1.05T C : pseudoscalar (upper panel) and vector (lower 
panel) channels. 
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FIG. 2: As in Fig. [T] but for the P-wave correlators: scalar 
(upper panel) and axial-vector (lower panel) channels. 



In Figs. [T] and [5] we display our results for 
G H (T,T)/G% c {T : T,T r ) for charmonium at T r = 1.05T C 
and Nf — 2. They should be compared to the lattice 
calculations of Ref. [|[ (their Figs. 3 and 4, respectively). 
They should also be compared to the results of Ref. [l5| , 
where the present model has been employed without in- 
cluding the zero mode contributions (note that in a pure 
potential model calculation the ratio G/G ICC depends 
only on the angular momentum state, S or P-wave, and 
not on the specific channel: thus, for instance, in Ref. fl5j 
the ratio in the vector channel coincides with the pseu- 
doscalar one). 

The temperature dependence of the ratio is now 
well reproduced in the scalar and axial- vector channels, 



thanks to the constant contribution of Eq. (|7b[) . In the 
vector channel the strong enhancement of the ratio at 
large r is also well described, owing again to the zero 
mode term in the spectral function. On the whole, the 
data of Ref. Q are described semi-quantitatively: we 
stress again, as thoroughly discussed in Ref. 15], that 
an accurate description of the lattice data should not be 
expected, due to the discretization effects in the lattice 
calculations. 

In the pseudoscalar channel, where no constant contri- 
bution to the correlators is present, the ratio evaluated 
on the lattice tends (also in the quenched case) to stay 
around one up to some temperature (interpreted as the 
dissociation one) and then it starts decreasing; in our cal- 
culation, on the other hand, the ratio is monotonically 
decreasing with r at all the temperatures. A remnant of 
this behavior is also visible in the vector channel, where 
the ratio drops slightly below one, before jumping up. 

Indeed, the binding due to the potential extracted from 
the lattice free energies gets softer with increasing T and 
this is reflected into the decrease of the square wave func- 
tion at the origin |l5|. On the other hand, the lattice 
spectral functions generally display ground state peaks 
of nearly constant strength up to the dissociation tem- 
perature. Note that when in the potential model the 
bound state energy is going to zero, a strong resonance- 
like pjlj contribution appears in the continuum (see the 
spectral functions displayed in Ref. 15]). Although it 
might be possible that on the lattice this contribution 
has been mistakenly identified as a bound state, in our 
calculations it still has a strength gradually decreasing 
with T. 

Note, anyway, that the trend predicted in the potential 
model is qualitatively correct, since the quenching of the 
pseudoscalar ratio gets stronger with increasing T. Simi- 
lar results are also obtained for the quenched calculations 
of Refs. 0, 0] and we do not display them here. 

In Ref. [HI an alternative explanation for the flatness, 
up to the dissociation temperature, of the G/G rcc ratio 
in the pseudoscalar channel has been proposed: in that 
potential model no continuum contribution is present 
when there are S'-wave bound states; thus, assuming 
the existence of just one bound state, one should have 
Gp S (r,T)/G^ s (r,T,T r ) = \R PS (0, T)\ 2 /\R PS (0, T r )\ 2 ■ 
K(T,M(T),T)/K(T,M(T r ),T) (note that this expres- 
sion in general is not normalized to one at r = 0). How- 
ever, in Ref. [131 ] the ratio between the wave functions 
at the origin has been neglected: it has a strong tem- 
perature dependence and it is precisely this term that 
generates the quenching with increasing T of the ratio in 
the pseudoscalar channel. Since the mass of the bound 
states turns out to be quite stable, the remaining fac- 
tor, K(T,M(T),T)/K(T,M(T r ),T), is nearly constant 
and yields the flat behavior of the lattice calculations, 
which thus here appears just as an artifact of the ap- 
proximation. 

A potential model calculation including the zero modes 
has also been done in Ref. [lq . using a hybrid model 
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where the spectral functions are calculated by solving a 
Schrodinger equation with a screened potential at low 
energies and employing the perturbative expressions at 
large energies. In spite of the different models, also this 
approach yields very good results for the ratio G/Grec 
of the P-wave correlators, indicating that this ratio is 
essentially dominated by the zero mode contribution. On 
the other hand, the results of Ref. [l6| in the pseudoscalar 
channel are rather different from ours: in the cc system at 
T slightly above T c , the authors of Ref. [l6| get a rather 
flat ratio, compatible with the lattice data; when T grows 
the lattice ratios get increasingly quenched, while their 
calculations show a moderate enhancement (cf. Fig. Q] 
here and Fig. 11 of Ref. [HI]). According to the authors of 
Ref. [l(| the discrepancy should be acribed to relativistic 
effects. 

It should be noted that in Ref. [16] the finite T corre- 
lators are compared to the T — one and in the latter an 
ad hoc factor is introduced to account for radiative cor- 
rections: in the pseudoscalar channel this factor yields a 
100% correction to the low energy part of the spectral 
function. Apparently, no radiative corrections have been 
included at finite T. 

To summarize, we have updated a recent calculation 
[If| of quarkonium Euclidean correlators at finite T, by 



including the zero mode contributions that are present 
in the lattice spectral functions. The model developed 
in Refs. [H, [l|f is based upon an effective potential ex- 
tracted from QQ free and internal energies measured on 
the lattice. Good agreement with the lattice results has 
been found for the P-wave correlators (scalar and axial- 
vector channels). Also the vector channel appears to be 
well described. 

The pseudoscalar correlators, on the other hand, show 
the same T-dependence as the ones calculated on the 
lattice, but not the same T-dependence below the disso- 
ciation temperature. This discrepancy has already been 
discussed in Ref. [l5| as being related to the fact that the 
strength of the bound states (and of the zero-energy res- 
onances) in the potential model is gradually decreasing 
with T, whereas the peak associated to the ground state 
in the lattice measurements seems to keep a constant 
strength up to dissociation. Hopefully, more accurate 
lattice calculations should shed some ligth on this issue 
in the near future. 

We would like to stress that these results have been ob- 
tained without any adjustment of any parameter, the ef- 
fective potential being fixed by independent lattice mea- 
surements of the QQ free energies. 
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Although one in principle would like to use the physi- 
cal quark masses, it should be remembered that the lat- 
tice free energies are calculated for static (infinite mass) 
heavy quarks. 

Strictly speaking there are no S-wave resonances for po- 
tentials without a repulsive barrier; however, the large 
enhancements in the cross-section associated to ne arly 
bound states are often called zero-energy resonances [29| ] . 



